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ABSTRACT 

The double or extended main-sequence turnoffs (MSTOs) in the color-magnitude diagram 
(CMD) of intermediate- age massive star clusters in the Large Magellanic Cloud are always in- 
terpreted as age spreads. However, it has been confirmed that the age spreads do not exist in 
young massive star clusters in the Large Magellanic Cloud, which is challenging the scenario of 
the age spreads. The result of effects of rotation on the MSTOs of star clusters is con- 
flicting in previous works. Compared with previous works, we consider the effects of 
rotation on the MS lifetime and take a different efficiency of rotational mixing. Our 
calculations show that rotating models have a fainter and redder turnoff with respect to non- 
rotating counterparts with ages between about 0.8 and 2.2 Gyr, but they have a brighter 
and bluer turnoff when age is larger than 2.4 Gyr. The spread of the MSTO caused by a typical 
rotation rate is equivalent to the effect of an age spread of about 200 Myr in the CMD 
of intermediate-age star clusters. Rotation could lead to the double or extended MSTOs in the 
CMD of the star clusters with ages between about 0.8 and 2.2 Gyr. But the extension is 
not significant, even does not exist in young clusters. If the efficiency of the mixing were 
high enough, the effects of the mixing would counteract the effect of the centrifugal 
support in the late stage of evolutions, the rotationally induced extension would 
disappear in the old intermediate-age star clusters; and young clusters would have 
a blueward extended MSTO. 

Subject headings: stars: rotation — stars: evolution — globular clusters: general 



Introduction 



Since JMackev fc Brobv Nielsen] (|2QQ7h and 
Mackev et al.l (J2008I ) discovered the double 
main-sequence turnoffs (MSTOs) in the color- 
magnitude diagram (CMD) of intermediate-age 
star clusters, such as NGC 1846, 1806, and 
1783 in the Large Magellanic Cloud (LMC), the 
phenomenon of the double or extended MSTOs 
have been dis covered in th e mor e and more 
star c lusters dMackev et all l2008t JGlatt et al.1 



Goudfrooii et all 120091 l2011aL iKeller et all 12012 : 
Piattil |2013[ ). Some star clusters, s uch as NGC 



1751 in the LMC and NGC 419 (iMiloneejaL 



20091 : IGirardi et all 120091 : iRubele et all 1201 ll ) 



20081 : iMilone et al.l 120091 : IGirardi et al.l 2009; 



m 

the Small Magellanic Cloud, have not only the 
double or extended MSTO s but two distinct red 
clumps (RCs). Recently. [Girardi et"aD <|2013h 
found the NGC 411 in the Small Magel- 
lanic Cloud also has an extended MSTO. 
Moreover, Omega Centauri ex hibits a main 
(MS) bifurc ation dPiotto et al 
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and NGC 2808 dPiotto et all E 
and NGC 6752 (IMilone et alJ l2013h pos- 
sess a triple MS split. These discoveries are 
not consistent with the classical hypothesis that 
a star cluster is composed of stars belonging to 
a simple, single stellar population with a uniform 



age and chemical composition. 

The double or extended MSTOs have 
mostly been interpreted as that the population of 
those star clusters have bi m odal a g e distributions 
dMackev fc Brobv Nielsen! l2007t iMackev et aL 
200 81) or age spreads of about 100 - 500 Myr 



MS will be prolonged, vice versa ( Maeder 



■ Mil one et aL 2009; G irardi et al. 2009; Ru bele et al 
2Q10L [201 it iKeller et al.ll2012l : IPiatti 2013). How- 



Platais et al.l (J2012I ) noted that this long 



ever, 

period of star formation seems to be odds with 
the fact that none of the younger clust ers are 
known to have such a trait. iBastian fc Silva-Villal 
(|2013l ) found that young massive clusters NGC 
1856 and 1866 in the LMC do not have such 
large age spreads, which strongly queries the 
scenario of the age spreads. The age is 



abou t 280 Myr for NGC 1856 (jKerber et al 

160 



2007h and ~ 160 - 25 Myr for NGC 1866 

([Brocato et alJ 120031 ). Moreover, the lasting 



time for explaining the double MSTOs and the 
dual R Cs of the NGC 17 51 seem to be not con- 
sistent (JYang et al.l 120 111 ). The contributions of 
interactive binary stars to t he double MSTO s 
and dual RCs were studied bv lYang et al.l (|201l[ ). 
They found that binary interactions and merging 
can reproduce the dual RCs and extended MSTO 
in the CMD of an intermediate-age star cluster. 
However, the fraction of the interactive binary 
systems is too low to explain the observed prop- 
erties. Moreover, the binary interactions make 
the stars have a younger apparent age than non- 
interactive systems, i.e. young stars are minor- 
ity, which seems to be contrary to the finding of 
Milone et aD ([20091 ). who found that about 70% of 



stars belong to the blue and bright (young) MSTO 
and around 30% to the red and faint (old) MSTO. 
The predict of the interactive binaries seems to be 
contrary to this finding. 

The effects of rotation on the structure and evo- 
lution of MS stars are mainly in three ways: (1) 
The effect of the centrifugal support leads to 
a decrease in the effective temperature and lumi- 
nosity by decreasing the effective gravity; (2) the 
mixing of elements results in a decrease in the stel- 
lar radius and an increase in the effective tempera- 
ture by changing the distributions of elements; (3) 
the Von Zeipel effect and the mixing of elements 
can affect the instability of convection by changing 
the radiative and adiabatic temperature gradients. 
If the convective core increases, the lifetime of the 



Talon et al 



1997; Maeder & Mevnetl 



19971; 
2013ah 



1987; 



2000; 



Yan g et al 

The effects of rotation on the MSTO of 
interm ediate- age s tar clu sters were first stud- 
ied bv IBastian fc de Minkl (|2009h . They obtained 
that rotating stars have a lower effective tem- 
perature and luminosity than non-rotating coun- 
terparts and concluded that stellar rotation can 
mimic the effect of a double population. However, 
iGirardi et al.l ()201ll ) calculate d the evolution of 
rotating models by using the Eggenberger et al. 
(|2010f ) code and obtained that rotating models 
have a slightly hotter and brighter turnoff with re- 
spect to non-rotating ones. They concluded that 
rotational effect could not e xplain the presence 
of the extended MSTO. In IBastian fc de Minkl 
(|2009f ) models, the effect of rotation on the age of 
stars was neglected and the effects of the rotation- 
ally induced mixing may be w e aker compared 



to those in IGirardi et al.l ( 201ll ) models 



which leads to the fact that their rotating mod- 
els have a lower temperature and luminosity than 
non-rotating ones in the whole stage of the MS. 
Moreover, they only calculated the evolution of 
the star with a mass of 1.5 Mq, assuming that 
the results for the star can be applied to others. 
In fact, the effects of rotation on the structure and 
evolution of stars depe nds on the mass of stars for 
a given rotation rate (IMaeder fc Mevnetl l2000l ; 
Yang et al.l l2013ah. On th e contrary, in the 



models of iGirardi et al.l (j201ll L the effects of the 
mixing may be more efficient, which results in 
the fact that their rotating models have a higher 
effective temperature than non-rotating ones in 
the late stage of the MS for stars of any mass. 



In the early stage of the MS of lEggenberger et al. 



(2010) models, the rotating models also exhibit a 
lower effective temperature and luminosity, which 
mainly results from the effect of the centrifugal 
support. The latest evolu t ionary tracks cal- 
culated by Georgy et al.l (J2013I ) also show 
that the rotating models have a lower ef- 
fective temperature compared to the non- 
rotating ones in the early stage of the MS 
(see their Fig. 11). I f the extent of the eff ect of 
mixing in the models of Gi rardi et al.l (j201lf ) were 
different, the results might be changed. Thus the 
effects of rotation on the MSTO need to be care- 
fully rechecked. 



Furthermore, the observed characteris- 
tics of uj Centauri and NGC 2808 can be in- 
terpreted by that a fraction of globular clus- 
ter stars have sizable helium e nhancements 
over the primordial value ( iPiotto et al 
2005L l2007h . In addition, there are abun- 



Yang fc Bil (|2QQ7h and lYang etaL 



The 



dance anomalies in o t her globular cluster s 



(JGratton et alJ 120041 : iMilone et all l2013h . 



Vent ura et al 
(J200 



S) 



In order to understand the self-enrichment 
in globular cluster stars, 
(|200lh and iD'Antona et al 
gest that the low-mass stars 



sug- 
have 



may 

been polluted at the surface by accre- 
tion from the gas that was lost from the 
evolving intermediate-mass asymptotic gi- 
ant branch stars, which requires a timescale 
of a few 100 Myr. However, mass shed- 
ding models from fast ro tating massive 
stars (Decressi n et alJ 2007) o r massive bi- 
naries ( De Mink et al.l 1200 9) can enrich 
the cluster on a timescale of a few Myr. 
The rotational mixing in massive stars 
(|Hunter et al.ll2008l: iBrott et aDl201ll) and 



solar models (lYang fc Bil 120071 : iBi et al. 
20111 ) is far from settled. Other mixing 
processe s may exist in the Sun and stars 

20081 : 



Hunter et al 



(Basu & Antia 2008; 

IBrott et al.N201lh , such as magnetic fields, 
gravity waves, and so on. 

The evolutionary tracks of rotating mod 
els with different rotatio n rates were pub- 



lished by several groups (jBrott et al.l 2011; 



Georgy et al. 2013). However, their calcu- 
lations did not cover the mass range be- 
tween 1.3 and 1.7 M . In this work, we also 
focus on the effects of rotation on the MSTO of 
intermediate- age massive star clusters. The pa- 
per is organized as follow. We simply show our 
stellar evolutionary tracks in section 2. The re- 
sults are represented in section 3. In section 4, we 
discuss and summarize the results. 

2. Evolutionary tracks 

2.1. Evolution code and assumptions 

We used the Yal e Rot ation Evolution Code 

HIS 



( Pinsonneault et al. 



Yang fc Bil 12007 ) to 



compute the evolutions of stellar models with and 
without rotation. The input physics and some 
initial parameters are the same as those used in 



meridional circulation, the Goldreich-Schubert- 
Fricke instability, and th e secular shear instability 
are co nsidered (Endal fc Sofia! 19 78c IPinsonneault et al.l 



1989) in all rotating models. The magnetic effects 
(|Yang fc B i 2006) are only considered in stars with 
a convective envelope, i.e. with mass less than 
about 1.3 Mq. In more massive stars, the mag- 
netic effects, magnetic braking, and mass loss 
are neglected, i.e., the total angular mo- 
mentum is assumed to be conserved. The 
transport of angular momentum and the mixing 
of elements are treated as a diffusion process in 
our m odels ([Pinsonneault et al.l 1 1989c I Yang fc Bi 
20061 ). The efficiency of rotational mixing 



is described by a parameter / c , which is 
less than unity and used to account for 
how the instabilities mix material less effi- 
ciently than they transport angular momen- 
tum. Usually, the value of f c is about 0.02- 
0.03 riChabover fc Zahnl Il992l: lYang & bH 
20061 : iHunter et al.ll2008l : IBrott et al.ll201lh 7 
Here, we adopted a higher value, i.e., 0.2. 

The initial metallicity Z and hydrogen abun- 
dance X was fixed at 0.008 and 0.743, respec- 
tively. For a given mass, the rotating and non- 
rotating models share the same initial parame- 
ters except for the rotation rate. Rotating mod- 
els were evolved from the zero-age MS (ZAMS) to 
the end of the MS, assuming the initially uni- 
form rotation of 1.0 x 10 -4 , 1.5 x 10 -4 , and 
2.0 x 10 -4 radian s _1 which corresponds to 
an initial period of about 0.73, 0.49, and 
0.37 day or a rotation rate of about 0.2, 
0.3, and 0.4 times Keplerian rotation rate 
on ZAMS, respectively. The period of 0.49 
day produces a surface velocity of about 130-180 
km s _1 for stars with masses between 1.3-3.0 M 
on Z AMS, whose typica l velocity is about 150 km 



([Rover et al.l 12007). Compared with the 



initial velocity on the ZAMS of 150 km s _1 
in Girardi et al. (2011) models, our initial 
velocity i s varing with mass . In th e simula- 
tion of iBastian & de Minkl (|2009h . the ini- 
tial rotation rate is around 0.4 times Kep- 
lerian rotation rate on the ZAMS, which is 
higher than ours. The value of f c is 0.0228 
in iBastian fc de Minkl (2009) models (see 
IBrott et alJ (|201lh ). which is lower than our 
value of 0.2. 



2.2. Evolutionary tracks 

The evolutionary tracks of rotating and non- 
rotating models with M = 1.15 and 2.2 M are 
shown in Fig. [TJ For the star with M = 1.15 M , 
due to the presence of magnetic braking, the ef- 
fects of mixing are almost dominant in the whole 
MS stage. Thus the rotating model mainly ex- 
hibits a higher effective temperature than non- 
rotating one at the same age. For stars with 
mass less than 1.3 M , magnetic braking 
leads to the fact that rotation rate decreases 
rapidly in a few hundred Myr. Thus the 
effects of rotation on stars are dominated 
by the rotational mixing, even more sensi- 
tive to the efficiency of the mixing than to 
the rotation rate. For the star with M = 2.2 
M , in the early stage of the MS phase, the ef- 
fect of centrifugal support dominates the influ- 
ences of rotation on the structure and evolution of 
the star; therefore the rotating model has a lower 
effective temperature and luminosity than non- 
rotating one. However, as the evolution pro- 
ceeds, efficient mixing leads to the fact that 
the outer edge of convective core is slightly 
extended (see Fig . [2]), which enhances the 
effects of mixing dMaeder fc Mevnetl l2000l : 



Yang et al.l l2013al ) . The efficient mixing pro- 
longs the lifetime of the MS by feeding fresh hy- 
drogen fuel into the hydrogen-burning region and 
makes the rotating models exhibit a higher ef- 
fective temperature than non-rotating one at the 
same age, wh ich is consistent with the resu lt of 
Maederl (|l987h and lEggenberger et al.l (|2010h . 

The evolutionary tracks of models with M = 
1.4, 1.5, and 1.6 M are shown in Fig. [3j where 
the positions of the models with a given age are 
labeled by the notation i o\ In the early stage 
of the MS, due to the fast rotation and the fact 
that the gradients of chemical compositions are 
small, the effect of rotational mixing is not signifi- 
cant and the effect of the centrifugal support 
plays a dominant role, thus the rotating models 
have a lower luminosity and effective temperature 
than non-rotating one. But as the evolution pro- 
ceeds, the rotating models can exhibit a higher 
luminosity and a lower effective temperature than 
non-rotating one at the same age. The effective 
temperature of rotating models also can be higher 
than that of the non-rotating one, which depends 
on the age and rotation rate of stars. 



The effect of the centrifugal support leads 

to a decrease in the effective temperature and 
gravity of stars. The decrease of the gravity acts 
as an effective decrease in the mass of stars. 
The mass of the convective core of MS stars de- 
creases with decreasing stellar mass. Therefore, 
the centrifugal effect results in a decrease in the 
convecti ve core ( see the Fig. [2]). The simula- 
tion of Ijulienl ( 19961 ) also shows that con- 
vective penetration may be hindered in ro- 
tating stars. Hydrogen is ignited at the same 
temperature. The smaller the convective core, 
the less the hydrogen that can be consumed in the 
core, the shorter the lifetime of the MS. Hence 
the effect of centrifugal support leads to an 
acceleration of the evolution of stars. For exam- 
ple, the rotating model with M — 1.5 M 
and P = 0.49 day has evolved to the MS 
hook at the age of 1.5 Gyr. However, the 
non-rotating counterpart is still on MS with 
a higher effective temperature and a lower 
luminosity (see Fig. [3]), which leads to that 
the rotating model seems to be older than 
non-rotating one. In fact, they have the 
same age. 

The rotational mixing can bring hydrogen fuel 
into the core from outer layers and transport the 
products of H-burning outwards. The helium in 
the convective core is brought to the radiative re- 
gion, which leads to an increase in the density p 
at the bottom of the radiative envelope. The adi- 
abatic gradient V a d is proportional to 1/p, thus 
the convective core can be slightly increased by 
the mixing. The increase bring more hydrogen 
fuel into the core and enhance the effect of mixing 
which prolong the lifetime of the MS of stars and 
lead to an increase in the effective temperature by 
increasing the mean density of stars. 

The effects of rotation on the structure and evo- 
lution of stars are a result of the competition be- 
tween the effect of the centrifugal support 
and that of rotational mixing. For stars with 
mass larger than 1.3 M , due to without magnetic 
braking, the effect of the centrifugal support 
plays a dominant role in the early stage of the 
MS. In the late stage of the MS, although the ef- 
fect of the centrifugal support plays an impor- 
tant role, the effects of mixing partly counteract 
the influences of the centrifugal effect, even are 
dominant. Thus, the convective core of the ro- 



tating 1.5 M model with the initial rotation 
period of 0.49 day is almost as larger as that of 
non-rotating one in the late stage of the MS (Fig. 
[2]). In the case of the evolution of rotating 
model with M = 1.5 M and P = 0.37 day, 
the effects of mixing counteract the effect 
of the centrifugal support in the late stage 
of the MS. Thus the luminosity and effec- 
tive temperature of the rotating model are 
approximate to those of non-rotating one at 
the age of 1.5 Gyr. But before the age of 
about 1.2 Gyr, the rotating model has an 
obviously lower effective temperature com- 
pared to the non-rotating one. Our calcu- 
lations show that all stars with Z = 0.008 and 
masses between around 1.3 and 2.0 M have 
similar characteristics. 

We also calculated the evolutionary tracks 
of rotating models with an initial rotation 
period of 0.49 day and a low efficiency of 
mixing (f c = 0.03). Figure [4] shows the evo- 
lutionary tracks of models with M = 1.15 
and 1.5 M . For a comparison, the track of 
the rotating model without mixing (f c — 0) 
of elements but with other effects of rota- 
tion, such as the effect of the centrifugal 
acceleration, is also plotted in the panel of 
the star with M = 1.5 M . The lower the 
efficiency of mixing, the lower the effects 
of the mixing. Thus the rotating models 
with a small f c have lower effective tem- 
peratures compared to the rotating models 
with a large f c . 

3. Isochrones and synthesis results 

The rotating models of stars with masses 
between about 1.3 and 2.0 M have a lower 
effective temperature than non-rotating 
ones, which would lead to a spread in color of 
the intermediate- age massive star clusters. Thus 
we calculated a grid of evolutionary trackqj of 
rotating and non-rotating models with masses be- 
tween 1.0 and 3.0 M . The mass interval SM 
is between 0.01 and 0.02 M for stars with 
masses between 1.15 and 2.0 M but is about 
0.1 M for others. The metallicity Z of evolu- 
tionary models was converted firstly into [Fe/H]. 
Then the theoretical properties ([Fe/H], T e ff, logg, 



log L) have been transformed into colors and mag- 
nitudes using the c olor transformation tables of 
Leieune et all (| 19981 ). 



lr rhe tracks can be obtained by e-mail to Wuming Yang. 



Figure [5] shows the CMDs of different isochrones 
obtained from our evolutionary models. The all 
rotating models have the same initial rota- 
tion period of 0.49 day. The MSTO of rotating 
models with age = 1.2 and 1.5 Gyr is almost coin- 
cident with that of non-rotating models with age 
= 1.35 and 1.7 Gyr, respectively. This indicates 
that the rotationally induced spread in color is 
equivalent to the effect of an age spread of 
about 200 Myr. With the increase or de- 
crease in age, the extension becomes nar- 
rower and narrower. For example, when 
the age increases to 1.7 or decreases to 0.9 
Gyr, the extension is equivalent to the effect 
of an age spread of about 100 Myr. This 
is because the effect of the centrifugal ac- 
celeration is partly counteracted by the ef- 
fects of mixing. When the age is located 
between 0.8 and 1.9 Gyr, rotating models ex- 
hibit a redder and fainter turnoff with respect to 
their non-rotating counterparts, which leads to 
the fact that rotating population are appar- 
ently 100-200 Myr older than non-rotating 
counterparts in the CMDs. However, they 
have the same age. The evolutions of rotating 
models with masses between about 1.3 and 2.0 M 
are faster than those of non-rotating ones. Thus 
the mass of the turnoff stars of the isochrone of 
rotating models is slightly less than that of non- 
rotating ones. Therefore, the MSTO of rotating 
models is fainter than that of non-rotating ones. 

When the age of star cluster is older than 
2.4 Gyr, the rotating models exhibit a bluer and 
brighter MSTO with respect to their non-rotating 
counterparts; the rotationally induced spread 
is equivalent to the effect of an age spread of 
about 200-400 Myr. When the age of star 
clusters is younger than 0.6 Gyr, the rotat- 
ing models also exhibit a bluer and brighter 
MSTO with respect to their non-rotating 
counterparts. For example, in a star cluster 
with an age of 0.4 Gyr, the spread caused 
by the rotation is simil ar to the effect of an 
age spread of 50 Myr. Basti an &; Silva-Villal 
(J2013 ) found that the age spread is less than 
35 Myr in NGC 1856 and NGC 1866. When 
the mass of stars is larger than about 2 



Mq or less than about 1.3 M , the rotat- 
ing models mainly exhibit a higher effective 
temperature than non-rotating ones at the 
same age. Thus the rotating models have 
a bluer and brighter turnoff compared to 
their non-rotating counterparts when age is 
younger than 0.6 Gyr or older than 2.4 Gyr. 

In order to understand whether rotation can 
produce the extended or double MSTOs or not, we 
performed a stellar population synthesis by Monte 
Carlo simulation, for a total mass of about 3 x 10 4 
Mq (the mass range of the intermediate-age 
massive star clusters is between abou t (1- 
20)xl0 4 M© (JGoudfrooij et al.l IJOllbl )) fol- 
lowing the log -normal initial mass function of 
IChabrier (120011). I n orde r to match the result 
of iMilone et al.l (J2009h that about 70% of 
stars belong to the blue and bright MSTO 
and around 30% to the red and faint MSTO, 
we simply assumed that 70% of stars do not 
rotate but 30% of stars rotate with the same 
initial period of 0.49 day. This assumption 
may do not represent the realistic situ- 
ation. In the synthesized population, we in- 
clude observational errors, taken to be a Gaus- 
sian distribution with a standard deviation of 0.01 
an d 0.015 in magnitude and color as that 
of iBastian fc de Minkl (J2009h . respectively. 
The deviation of 0.015 in V-I corresponds 
to a deviation of about 90 K in the effective 
temperature in our models. 

The CMDs of the synthesized star clusters with 
different ages are shown in Fig. [6] The double 
MSTOs can be clearly seen in the star clusters 
with age = 1.2 and 1.5 Gyr, in which the hot and 
bright MSTO is dominant. The continuously ex- 
tended MSTO exists in the star cluster with age 
= 0.9 Gyr, where the hot and bright MSTO is 
also dominant. Compared to the extension of the 
MSTO of star clusters with ages between 0.9 
and 1.6 Gyr, that of the star clusters with 
age = 0.4 and 1.7 Gyr is not obvious. More- 
over there is almost no extension of the MSTO of 
the star cluster with age = 2.0 Gyr. For the star 
clusters with age = 2.6 and 3.0 Gyr, their MSTO 
is slightly extended by the effects of rotation, but 
the cool and faint MSTO is dominant. The dif- 
ference in the V-I between the isochrone 
of rotating models and that of non-rotating 
ones is less than 0.01 for the star cluster 



with an age of 0.4 Gyr. Thus the extension 
of the MSTO in this cluster is dominated 
by the deviation of 0.015 in color. However, 
the difference is about 0.04 for the star clus- 
ter with an age of 1.2 Gyr. The spreads in 
star clusters with ages between 1.2 and 1.5 
Gyr are mainly due to rotation. 

The isochrones of rotating models with 
an initial period of about 0.73 day and the 
synthesized results are shown in Figs. [7] 
and [H respectively. For the star clusters 
with age less than 0.7 Gyr, their MSTOs 
are not extended by rotating stars. When 
the age of star clusters is located between 
0.9 and 2.2 Gyr, the rotating models have a 
redder and fainter MSTO compared to the 
non-rotating counterparts, which is equiva- 
lent to the effect of an age spread of about 
100-200 Myr. Compared to the results ob- 
tained from the rotating models with the 
initial rotation period of 0.49 day, the re- 
sults for the star clusters with age larger 
than 2.6 Gyr are the same; but the red and 
faint MSTO appears in older intermediate- 
age star clusters. 

Figures [9] and [10] show the results ob- 
tained from the rotating models with an 
initial period of 0.37 day. The interesting 
scenarios are that the spread caused by ro- 
tation disappears in star clusters with ages 
between 1.4 and 2.0 Gyr, but the spread 
that is similar to the effect of an age spread 
of about 100-200 Myr still exists in the star 
clusters with ages between 0.8 and 1.3 Gyr. 
When the age is less than 0.7 Gyr, the rotat- 
ing models have a bluer and brighter MSTO 
compared to non-rotating ones, which is 
equivalent to the effect of an age spread of 
about 100 Myr. This is because the effects 
of mixing are easier to play a dominant role 
in the late stage of the MS of stars for a 
high rotation r ate. These re sults are simi- 
lar to those of iGirardi et al.l (|201ll ) except 
for star clusters with ages between 0.9 and 
1.2 Gyr. The MSTO of star clusters with 
age larger than 2.6 Gyr is similar to that 
of the star clusters with a lower rotation 
rate. The efficiency of mixing in our mod- 



els should be less than that of Girardi et al. 



(|201ll ). Thus, in our calculations, it needs a 



higher rotation rate to produce the similar 
results. 

The NGC 1751 has the double MSTO 



and i ts age is about 1300-1500 Myr (jMilone et al.l 
2009). In a star cluster, the initial rotation 
period of stars should be different. For 
simplicity, we assumed that its population 
is composed of 70% non-rotating stars and 
30 % rotating stars which is obtained by 
interpolating between our rotating models 
with different rotation rates, assuming the 
initial rotation rate has a Gaussian distri- 
bution with a peak at about 0.3 times Ke- 
plerian rotation rate on ZAMS and a stan- 
dard deviation of 0.06. This assumption 
may do not represent the realistic distribu- 
tion. For example, the distribution of rota- 
tion rate of A0-A1 type stars are more flat 
( Rover et al.l 120071 ). A distance modulus of 
18.5 is adopted in this simulation. The sim- 
ulated results are shown in Fig. 111^ where 
the double MSTO can be seen. The struc- 
ture of the simulate d MSTO is very si milar 
to the observed bv lMilone et al.l (J2009h (see 
their Fig. 9). 

Comparing the isochrones obtained from 
rotating models with different rotation 
rates, one can find that the largest exten- 
sion is almost the same, i.e., can be equiva- 
lent to the effect of an age spread of about 
200 Myr. With the increase in rotation 
rate, for the intermediate-age star clusters, 
the extension of the MSTO will disappear 
in the star clusters with an older age due 
to the fact that the effects of mixing coun- 
teract the effect of the centrifugal support; 
for the old-age star clusters, the blue and 
bright extension of the MSTO will be not 
affected because magnetic braking domi- 
nates the rotation of stars with mass less 
than about 1.3 M ; for young-age star clus- 
ters, an blue and bright extension of the 
MSTO can be produced by high rotation. 

The observed extension of the MSTO of 
the intermediate-age star clusters can be 
equivalent to the effect of an age spread of 
about 100-500 Myr, which is broader than 
the extension caused by rotating models 
with f c = 0.2. Figure [4] shows that the ef- 
fects of rotation on the evolution of stars 



are sensitive to the value of f c . Thus we 
calculated the rotating population with the 
initial period of 0.49 day but with a lower 
efficiency of mixing (f c = 0.03). The results 
for this calculation are shown in Figs. [T2l 
and 1131 When the age of star clusters is lo- 
cated between 1.4 and 1.7 Gyr, the exten- 
sion of MSTO caused by rotation is equiva- 
lent to the effect of an age spread of about 
400 Myr, which decreases with increasing 
or decreasing age. But the CMD of star 
clusters with age less than 0.6 Gyr is al- 
most not affected by rotation. When the 
age is larger than 2.6 Gyr, rotating models 
exhibit a bluer and brighter MSTO com- 
pared to non-rotating counterparts, which 
can be equivalent to the effect of an age 
spread of about 350 Myr in the cluster with 
an age of 3.0 Gyr. 

4. Discussion and Summary 



According to the st udies ofjWolff et al. (2004), 
Zorec fc Roverl (|2012h . and lYang et al.l (J2Q13bh . 
the ZAMS models of A-type stars should be a 
Wolff ZAMS model which is required for the 
studies of the evolution of surface v elocity and 
the tr ansport of angular momentum (Yan g et al.l 
2013b). For simplicity in the computation, we 
adopted the ZAMS model with a uniform rota- 
tion. The rotating models computed from the 
Wolff ZAMS model also exhibit a lower effective 
temperature and a higher luminosity with respect 
to non-rotating ones in the late stage of the MS of 
stars with masses between 1.3 and 2.0 M , which 
is consistent with the results obtained from the 
evolution of the ZAMS model with a uniform ro- 
tation (see Fig. HJ. Thus our results about the 
effects of rotation on the CMDs cannot be signif- 
icantly changed by the ZAMS model. 

For a given rotation rate, the coefficient of the 
diffusion increases with increasing mass, i.e. the 
efficiency of the mixing increases with increasing 
mass, which leads to the fact that there is a 
critical mass M c for the effects of rotation 
on the structure and evolution of stars. For 
stars with mass larger than M c , the effect of 
the mixing soon exceeds the centrifugal effect in 
the early stage of the MS. For stars with masses 
between about 1.3 M and M c , the effect of 



the centrifugal acceleration is dominant for a 
long time in the early stage of the MS; the more 
massive the star, the larger the fractional lifetime 
of the MS that is dominated by the effect of 
the centrifugal acceleration because the cen- 
trifugal effect also increase with increasing mass 
for a given rotation rate. Our ZAMS models 
with mass larger than 2.0 M have had a 
little gradient of elements in interior com- 
pared to lower mass stars in which the gra- 
dient can be neglected. In addition, the MS 
lifetime is longer for the lower mass stars. 
Thus the diffusion of elements takes place 
at the beginning of evolution for stars with 
mass larger than 2.0 M , but it acts in a 
later time for stars with mass less than 2.0 
M . Thus the early evolution of the MS of 
stars with mass less than 2.0 M is easier 
dominated by the effect of centrifugal ac- 
celeration. The critical mass is about 2.0 
M in our models, which is c lose to the 
critica l mass of the He- flash of I Yang et al. 



( 20121 ). It should be affected by the mixing 
processes. Thus the value of about 2.0 M 
is debatable. 

In rotating stars, meridional circulation 
is an advection process. The rotational 
mixing would be much more complex than 
that described by a diffusion process. The 
value of 0.2 of f c in our models is obvi- 
ously higher than that calibrated against 
the nitrogen abundance in massive stars 
(0.02 28) ( Hunter et~aT] 120081 ; iBrott et aD 
201 ll) or against the solar model (0.03) 
(JYang fc BiM2006h . There may be other 
mixing mechanisms undescribe d by present 
model in stars and Sun ([Hunter et al. 



20081 : iBasu fc Antial 120081 ). such as grav- 
ity waves which are effici ent at transport- 
ing angular momentum (IZahn et al.l Il997l : 
Talon fc Charbonnell 12005^ The value of 
0.2 for the f c corresponds to a high efficient 
mixing in stars, which results in a higher 
mean density for stars with mass less than 
2.0 M when they approach the end of MS, 
i.e., the models with a high efficient mix- 
ing have a lower radius, which aids in ex- 
plaining the evolution of surfa ce velocity of 
stars with mass less 2.0 M o (jZorec fc Royerl 
20121 : lYang et al.ll2013bl ) besides producing 



the extension of the MSTO that is equiv- 
alent to the effect of an age spread of 200 
Myr. However, if the mixing is too effi- 
cient, which will bring an amount of hydro- 
gen into the burning region of stars. As a 
consequence, the extension of the MSTO 
caused by rotation would disappear. Our 
current understanding of the mixing pro- 
cesses should be limited. And the evolu- 
tions of rotating stars should exceed the 
predictive power of one-dimensional rotat- 
ing stellar models. The mixing in stars is 
debatable. 

The rotationally induced extension of 
MSTO do not exist in the simulated clusters 
with age between 1.4 and 2.0 Gyr when the 
initial rotation period is 0.37 day, which 
is similar to the results of iGirardi et al 



(|201ll ). In addition, the extension also dis- 
appears in the young clusters with the ro- 
tating models with Pq = 0.73 day. If we 
took a fixed rotation velocity Vq as the ini- 
tial rotation parameter for all models, the 
high mass stars would have a lower rota- 
tion rate but the low mass stars would 
have a higher rotation rate. As a conse- 
quence, one could obtain the result that 
rotation cannot affect obviously the MSTO 
of star clusters. Moreover, the efficiency 
of mixing s hould be different between our 
models and lGirardi et al.1 (|201ll ) models be- 
cause there is no the adjustable param- 
eter fc in the Geneva stella r evolution 



code (JMaeder fc Mevnetl l200oh . The dif- 
ferent initial rotation parameter and efR 
ciency of mixing would lead to the d iffer 



ences between our results and Girardi et al. 



(|201ll Vs. 

Compared with lBastian fc de Minkl (J2009h 
models, we calculated a dense grid of evolu- 
tionary tracks of rotating models and con- 
sidered the effects of rotation on the life- 
time of the MS in our isochrones and sim- 
ulations. Our results are similar to that of 
iBastian fc de Minkl (J2009h . i.e., the MSTO 
of star clusters with ages between about 
0.8 and 2.2 Gyr can be extended by stellar 
rotation. The extension is dependent on 
rotation rate, the efficiency of mixing, and 
the age of star clusters. 



For stars with masses between about 1.3 
and 2.0 M©, the effect of the centrifugal 
support plays a dominant role in the early 
stage of evolutions, which leads to the fact 
that rotating models have a lower effective 
temperature and evolve faster than non- 
rotating ones. As the evolutions proceed, 
the effect of rotational mixing partly coun- 
teracts the influence of the centrifugal ac- 
celeration, even dominates the effects of ro- 
tation on the structure and evolution of 
stars. As a consequence, the MSTO of 
intermediate-age star clusters with ages be- 
tween 0.8 and 2.2 Gyr can be extended red- 
ward, but the extension decreases with in- 
creasing or decreasing age. In the rotating 
stars with mass less than 1.3 M , the effect 
of rotational mixing plays a dominant role 
during the MS stage, which results in the 
fact that rotating models have a higher ef- 
fective temperature than non-rotating ones 
at the same age. Thus the MSTO of the 
older star clusters is extended blueward. 
The evolutions of rotating stars with mass 
larger than 2.0 M© are mainly affected by 
the effect of rotational mixing which de- 
pends on rotation rate and the efficiency of 
mixing. Therefore, the MSTO of young star 
clusters can be slightly extended blueward 
or cannot be extended, which is dependent 
on the rotation rate and the efficiency of the 
mixing. 

In t he 16 star clusters studied by lMilone et al. 
(J2009h (see their Table 3), the star clusters 
with ages between 950 and 1400 Myr have 
an extended MSTO of about 150-250 Myr, 
but the star clusters with age larger than 
1550 Myr have a smaller extension (Aage 
< 100 Myr), which seems to be similar to 
our results that rotationally induced ex- 
tension decreases with the increase or de- 
crease in the age of star clusters and to be 
more similar to the results obtained from 
the high-efficiency models than to those ob- 
tained from the low-efficiency models. 

The extension of the MSTO caused by 
rotation is significantly affected by the ef- 
ficiency of mixing and decreases with in- 
creasing the efficiency. The star cluster SL 
529 with an age of about 2.0 Gyr has an 



age spread of about 500 Myr (|Piattill2013h . 
The rotationally induced extension of the 
MSTO of star cluster with an age of 2.0 Gyr 
is similar to the effect of an age spread of 
about 250 Myr in our models with f c = 0.03 
(see Fig. I12[) . If we took a lower initial rota- 
tion rate, the spread that can be equivalent 
to the effect of an age spread of 400 Myr 
would appear in star clusters older than 1.7 
Gyr. Thus the spread of star cluster SL 529 
seems to be not conflicting with our results. 
In this work, we calculated a grid of evolution- 
ary tracks of rotating and non-rotating stars 
with Z = 0.008 and masses between 1.0 and 3.0 
M®. For a given rotation rate, the effects of rota- 
tion on the structure and evolution of stars are 
dependent on the mass and age of stars. For 
the initial rotation period of 0.49 day and 
the efficiency of mixing of 0.2, rotating 
models have a fainter and redder turnoff 
with respect to the non-rotating counter- 
parts when age is located between about 0.8 
and 1.9 Gyr. When the age is larger than 
2.4 Gyr, rotating models have a brighter 
and bluer turnoff with respect to the non- 
rotating ones. However, when the age is 
less than 0.6 Gyr, the rotation almost not 
affects the MSTO of star clusters. The ex- 
tension of the MSTO caused by the rota- 
tion can be equivalent to an age spread of 
about 200 Myr when the age is located be- 
tween about 1.2 and 1.6 Gyr. The extension 
decreases with the decrease or increase in 
age. When the initial rotation period de- 
creases to about 0.37 day, the extension of 
the MSTO caused by rotation disappears 
in star clustars with ages between about 
1.4 and 2.0 Gyr, young star clusters with 
age less than 0.7 Gyr have a blueward ex- 
tended MSTO; but star clusters with ages 
between 0.8 and 1.3 Gyr still have a red- 
ward extended MSTO. When the initial ro- 
tation period increases to about 0.73 day, 
the MSTO of star clusters with age less 
than 0.7 Gyr cannot be affected by rota- 
tion, but the MSTO of star clusters with 
ages between 0.9 and 2.2 Gyr is extended 
redward. When the value of the efficiency 
of mixing decreases to 0.03, the extension of 
the MSTO caused by rotation can be equiv- 



alent to the effect of an age spread of about 
400 Myr. The rotationally induced exten- 
sion of the MSTO of intermediate-age star 
clusters is dependent on the rotation rate, 
the age of star clusters, and the efficiency of 
the mixing of elements. But the blueward 
extension of the MSTO of star clusters with 
age larger than about 2.4 Gyr is mainly de- 
pendent on the efficiency of the mixing. 

Our simulation shows that rotation could 
lead to the double or extended MSTO in 
the star clusters with ages between 0.8 and 
2.2 Gyr or larger than 2.4 Gyr. The rota- 
tion seems to cannot result in the extension 
of the MSTO of star clusters with age less 
than 0.6 Gyr; at least, the extension of the 
MSTO of these star clusters is not signifi- 
cant, compared to that of intermediate-age 
star clusters. Considered the difference in 
the initial rotation rates and the efficiency 
of mixing, the redward extended MSTO 
may be easier to be observed in star clus- 
ters with ages between about 0.9 and 1.5 
Gyr. 
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Fig. 1. — Evolutionary tracks of rotating and non- 
rotating models in the Hertzsprung-Russell (HR) 
diagram. The value of Pi(i = 1,2,3) is the initial 
period on the ZAMS. The notation V shows a 
position of the models with a given age. 
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Fig. 3. — Same as Fig. [T] but for models with 
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Fig. 5. — The CMDs of the different isochrones. 
The solid (green) lines show the isochrones of ro- 
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Fig. 6. — The CMDs of the synthetic star clusters 
at different ages. The total mass is assumed to be 
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log-normal initial mass function. We assumed 30% 
of stars rotating with Pq = 0.49 day and f c = 0.2. 
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Fig. 9. — Same as Fig. [5] but the initial period of 
rotating models is about 0.37 day. 
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Fig. 10. — Same as Fig. [6] but the initial period of 
rotating population is 0.37 day. 
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Fig. 11. — The CMDs of a simulated star cluster 
with the age of NGC 1751. A distance modulus of 
18.5 is adopted. Top panel shows the cluster with- 
out rotating population; bottom panel represents 
the cluster with 30% rotating population obtained 
by interpolating between rotating models with dif- 
ferent rotation rate. 



17 



0.4 Gyr, No Rot. 
_ 0.4 Gyr, Rot. 
. 0.6 Gyr, No Rot. 




■ 0.8 Gyr, No Rot. 



0.8 Gyr, Rot. 
0.95 Gyr, No Rot. 



0.2 0.4 
V-l 



0.2 0.4 0.6 0.8 
V-l 




1 .2 Gyr, No Rot. 



1.2 Gyr, Rot. 
1.5 Gyr, No Rot. 




1.7 Gyr, No Rot. 
_ 1.7 Gyr, Rot. 
. 2.1 Gyr, No Rot. 



0.2 0.4 0.6 0.E 



Z=0.008 




~C'~-- -•-■:-:---.-._ 




(j> 2.0 Gyr, No Rot. 




^^k 2.0 Gyr, Rot. 

X 2.25 Gyr, No Rot. 







■ 2.4 Gyr, No Rot. 
_ 2.4 Gyr, Rot. 

. 2.2 Gyr, No Rot. 



2.6 Gyr, No Rot. 
_ 2.6 Gyr, Rot. 
. 2.4 Gyr, No Rot. 



|\ 3.0 Gyr, No Rot. 

\ 3.0 Gyr, Rot. 

\ 2.65 Gyr, No Rot. 



0.6 0.8 
V-l 



1 



1 



Fig. 12. — Same as Fig. [5] but the value of f c is 
0.03. 



E > ; 


Z = 0.008,age = 0.6 Gyr ! 


( 


] 0.1 0.2 0.3 
V-l 


e > : 


- Z = 0.008,age=1 .2 Gyr . 





2 0.3 0.4 0.5 
V-l 


E *n - 


^~Z~0.008, age=1.7 Gyr 1 




0.3 0.4 0.5 0.6 
V-l 


E "n - 


Z=0.008, age=2.4 Gyr 1 




Fig. 13. — Same as Fig. [6] but the value of f c is 
0.03. 
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